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We present a new application and implementation of the so-called Pound locking technique for the interro- 
gation of superconducting microresonators. We discuss how by comparing against stable frequency sources 
this technique can be used to characterize properties of resonators that can not be accessed using traditional 
methods. Specifically, by analyzing the noise spectra and the Allan deviation we obtain valuable information 
about the nature of the noise in superconducting planar resonators. This technique also greatly improves the 
read-out accuracy and measurement throughput compared to conventional methods. 



I. INTRODUCTION 

Resonators are ubiquitous devices used in one form or 
another in most areas of physics and in many important 
applications, from filters in mobile phones to experiments 
in cavity Quantum Electrodynamics (QED) 1 . Over the 
past few years advances in on-chip quantum information 
processing have meant that planar microwave resonators 
have taken on new roles: as quantum buses^, for gen- 
erating Fock states^ and as read-out elements for solid- 
state qubits^. Superconducing resonators are also used 
as kinetic-inductance detectors^. These applications in- 
volve measurements that need to be done quickly and ac- 
curately. It is natural to look at techniques already in use 
in the fields of optics and frequency metrology for inspira- 
tion on how to meet these new measurement challenges. 
A commonly used technique for reading-out resonators in 
precision frequency metrology is Pound-locking^, rou- 
tinely used in quantum optics^ and precision frequency 
metrology^ but surprisingly little known outside those 
fields. 

The intrinsic properties of a resonator are described 
by its centre frequency v>q and unloaded quality fac- 
tor Q u — vq/Av, where Av is the width of the reso- 
nance. This in turn is related to the attenuation rate 
a = ■nv{j/2Q. All these parameters are routinely mea- 
sured using a vector network analyzer ( VNA) . Resonator 
measurements using a VNA are in principle accurate as 
long as the resonance shape is known and fitting tech- 
niques can be used to obtain Q and z/o, but are also quite 
slow, only measures the average behaviour and do not 
provide any direct information about the level of noise 
in the system. Another common approach is to use ho- 
modyne detection techniques where the carrier is mixed 
down to DC (using a mixer )». Pound-locking has a few 
important advantages over both these methods. The 
main improvement comes from the fact that the mod- 
ulation frequency f m (see below) can be chosen to be 
much higher than the 1/f corner frequency of amplifiers 
etc; this is important since microwave amplifiers are very 
noisy at low offset frequencies, meaning the noise in a 
homodyne measurement^^ scheme can be easily domi- 
nated by system noise. We typically use f m of 1-5 MHz, 
well above the corner frequency of our amplifiers. It is 



also a " real time" technique (within the bandwidth of the 
loop) in that it always "tracks" the centre frequency of 
the resonator irrespective of the magnitude of the fluctua- 
tion; leading to a significant speedup in the measurement 
throughput. Finally, it also allows for accurate measure- 
ments of very small shifts in the centre frequency of a 
resonator which is important in many experiments. 

The purpose of this paper is to describe a new ap- 
plication and appropriate implementation of this well- 
established technique. We have adapted Pound-locking 
to the read-out and characterization of superconducting 
microresonators. We will also describe how one can ob- 
tain information about the nature of the noise in samples 
by using analysis techniques from precision frequency 
metrology. It is our hope that this work will expand 
the applications of the Pound locking technique. 



II. POUND-LOCKING 



Figure 1(a) show the basic schematic of a Pound- loop. 



Here we will describe how this is implemented in the case 
of microwave oscillators; for a very intuitive discussion 
of how Pound locking is used in optics see the review 
by Blackii. The main building blocks are a voltage- 
controlled oscillator (VCO) that generates a carrier of 
frequency / c , a phase modulator used to add sidebands 
±f m to the signal, the resonator, a square-law power de- 
tector, a down-converting mixer and finally a regulator 
circuit which controls the VCO. The principle of opera- 
tion is as follows: assuming the carrier is near resonance 
and the modulation frequency has been chosen so that 
/,„ Av the sidebands do not enter the resonator (note 
that this means that no extra power is injected into the 
resonator; nor is the bandwidth limited by a). The phase 
modulated signal going into the sample will be of the 
form 11 



V m = V [J (P) + 2JM0) sin27r/ m i]e 



2njf c t 



(1) 



where Jo,i are Bessel functions and j3 is the modulation 
depth, it can be shown that the latter should be chosen 
so that the power in the sidebands is -3 dB smaller than 
the carrier. The resonator can be anything that has a 
"notch" type frequency response, most Pound loops are 
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(a)The basic Pound loop 



(b)Practical implementation of the Pound loop. Various auxiliary 
components (filters etc.) have been left out for clarity. 



FIG. 1. 





FIG. 2. Micrograph and equivalent circuit for one of our 
lumped element resonators. The electromagnetic radiation is 
primarily inductively coupled to the inductive part of the res- 
onator. The coupling strength is set by the mutual inductance 
between the coplanar transmission line and the resonator. 



implemented using a circulator (meaning the sample is 
actually probed in reflection) but our planar resonator 
(see below and ref J^) have this response in transmission, 
see fig. HJ The scattering parameter S^i for a piece of 
transmission line shunted by a resonator can be written 



S21 — 



2(l + 2jQ u y) 
g + 2(l + 2jQ u y) 



(2) 



Where g is a parameter determining the coupling 
strength and y is the normalized centre frequency (/ — 
v o)l v o- After going through the sample the signal is 
detected by a square-law diode after which it is de- 
modulated using a mixer. The circuit therefore creates 
an error signal that is due to the interference between 
the transmitted signal and the energy leaking out from 
the resonator back into the transmission line. The signal 
from the diode can be writte n 13 i 14 



oc 2V 2 kJ (p) Ji(/3){Im[5 2 i(/ c )](Re[5 2 i(/ c + f m )] 



Re[S 21 (/ c -/ m )])- 

Re[5 2 i(/ c )](Im[S 2 i(/ c + f m )\ + Im[5 21 (/ C - /„)])} (3) 



where e is the error signal; and k is a constant that de- 
pends on the conversion efficiencies of the diode and the 
de-modulator. The net effect of the Pound loop is that 
whenever the carrier is off resonance, the transmitted sig- 
nal is phase shifted and will have an AM component, i.e. 
there is partial PM to AM conversion. The AM signal is 
then detected and used to correct the error. Hence, ide- 
ally this circuit will "lock" the frequency of the carrier 
f c to the resonance frequency . Figure [3] shows what 
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the error signal e will look like for a loop with a typical 
resonator. Once the system has locked to the resonator 
it will only undergo small excursions from v§ and e is 
nearly linear, using eqJ2]we can simplify eq|3] 



kV n 2 



(3 + 2) ; 



(4) 



where we have assumed optimum modulation depth (/3 = 
1.08) and coupling (g = 2). Hence, the slope of the error 
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FIG. 3. Simulated typical error signal. The parameters are 
f m = ^o/2000, Q=50000,g=2 and (3 = 1.08. 



signal depends on the uncoupled quality factor and the 
coupling strength of the resonator. 



III. EXPERIMENTAL 

The a ctual set up used in our experiment can be seen 
in 1(b) Our resonators are thin- film devices microfabri- 
cated on 5x10 mm 2 chips (see refi 2 " for a full descrip- 
tion). The chip is glued to a copper cold-finger and 



bonded to a connectorized chip carrier which is enclosed 
in an aluminium box; the box is then mounted on the 
mixing chamber of a dilution fridge which has a base 
temperature of 30 mK. The in-signal to the resonator 
goes via a filtered and heavily attenuated (50 dB) line. 
The out-signal passes through an isolator and then into 
a low-noise cryogenic amplifier with a noise temperature 
of about 4 K. 

The microwave signal (hereon referred to as the carrier) 
is generated by a voltage-controlled oscillator, in practice 
this is usually a microwave generator where the output 
frequency can be modulated using an external voltage 
(at a rate of for example 10 kHz/V). Another option is 
to use a local oscillator (LO) held at a fixed frequency 
and mix it with an intermediate frequency (IF) from a 
low-frequency VCO (in practice e.g. a function genera- 
tor that can be frequency modulated) using for example 
a single sideband modulator (SSB). The performance of 
the latter setup is in principle better since a low phase- 
noise LO can be used (the phase noise of the IF source is 
negligible) but makes the circuit more complicated (and 
care has to be taken to minimize the amplitude of the 
side-bands). Another advantage of this latter arrange- 
ment is that the frequency shift can be read out directly 
with high accuracy using a standard counter. The car- 
rier then passes through a phase modulator which adds 
sidebands ±/ m . The shape of the error signal depends 
on the phase of the modulating signal, so this needs to 
be adjusted to compensate for lead/lag caused by the ca- 
bles, filters etc. The signal from the diode then goes to 
an RF lock-in amplifier (referenced to f m ) which is used 
with a very short time constant (as to not limit the band- 
witdh of the loop), essentially being operated as a mixer 
with gain. The error signal in the loop is the in-phase 



signal from the lock-in amplifier, and is fed to the PID 
controller. The output voltage from the controller then 
controls the VCO (which in our case outputs the IF) in 
order to null the error; we generally use a deviation of 
5-100 kHz/V. The PID controller is configured so that it 
attempts to keep the error signal as close to zero as possi- 
ble while still maintaining a loop bandwidth much larger 
than the frequency range of any noise we want to study. 
Care has to be taken to avoid dependence of the observed 
noise level on the loop parameters, this can happen e.g. if 
the gain is too high (resulting in oscillations that distort 
the noise spectrum) or the bandwidth too low(resulting 
in a roll-off of the noise at higher frequencies) . 

Most of our measurement data is acquired using a stan- 
dard 16-bit data acquisition system connected to the out- 
put of the PID controller. Sampling rates from 10 Hz to 
about 50 kHz are used, depending on the measurement 
time which varies from a few seconds to tens of hours. We 
also monitor the output signal using an FFT analyzer. In 
order to get an independent measure of the frequency we 
monitor the frequency of the IF source using a separate 
counter. All microwave equipment, counters etc in our 
experiment use an external 10 MHz frequency reference 
derived from a hydrogen maser. Our current setup allows 
us to track any change to the resonance frequency of the 
resonator with a bandwidth of a few kHz (limited by the 
PID controller) and a "real-time" frequency resolution of 
a few parts in 10 8 (limited by the intrinsic noise of the 
resonator). 

Our experiment is essentially dual to that of the con- 
ventional setup used to improve the stability of an oscil- 
lator. In a conventional experiment a "bad" oscillator is 
locked to a "good" resonator, we are here doing the ex- 
act opposite. The frequency noise level of our resonator 
is much higher than that of the carrier; by observing 
how the circuit strives to null the error signal we ob- 
tain direct information about the resonator. Moreover, 
since we are operating at offset frequencies smaller than 
/l = vq/2Q{~ 75 kHz), the Leeson effect 9 will not affect 
the shape of the frequency spectrum^. 

Our superconducting planar lumped element (LE) res- 
onators are about 300x300 /jm 2 in size, fabricated from 
niobium on a sapphire substrate. Several resonators 
(usually 5) are fabricated on the same chip and are cou- 
pled to the same transmission line allowing for frequency 
multiplexing. The data shown here comes from a res- 
onator with vq of 7.5 GHz and a loaded Q=50 000, it is 
similar to the ones discussed in ref 12 . Here we only show 
a subset of our results, for a full discussion seei&. 

In order to verify that our data really represents prop- 
erties of the sample and are not measurement artifacts; 
we also employ a number of reference samples. Here we 
will therefore also show data from a common dielectric 
resonator (DRO) in the shape of a 'puck' about 7 mm 
in diameter mounted inside a metal box and measured 
at room temperature. In this experiment we locked to 
a mode with a resonance frequency of 5.12 GHz and a 
loaded Q=5500. 
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FIG. 4. Root power spectral density (in units of Hz/V Hz) 
for a 7.5 GHz lumped element resonator measured at 30mK. 
the solid lines have a slope of -1 and -0.5, respectivly and are 
meant as guides for the eye. 

IV. DATA ANALYSIS 

The most straightforward way to measure noise is to 
use a spectrum analyzer or -equivalently- to calculate the 
power spectral density (PSD) from sampled time-domain 
data. Figure 0] shows a typical PSD from one of our LE 
resonators. Note that wha t is plotted is the frequency 
noise (in units of Hz/yHz) meaning the PSD represent 
spectral information of how the centre frequency of the 
resonator changes. The spectrum shows that random 
walk frequency noise (slope -1) dominates at low fre- 
quencies but that processes with a white (slope 0) spec- 
tra gradually takes over as the frequency increases. The 
solid line with slope -0.5 can be identified with flicker 
frequency noise. 

One limitation of PSD is that it is often very difficult 
to analyze "slow" processes. Another little known, but 
highly useful tool from the arsenal of frequency metrology 
is the Allan deviation (ADEV) 17 i 18 . Here we have used 
the overlapping Allan deviation a: 

^ M-2m+i (j+m-i "| 2 

° 2 = 2m 2 (M -2m+l) ^ [ ^ ^ l+m ~ ^ j 

(5) 

where M is the number of samples, To is the measurement 
timebase and y is the time-averaged fractional frequency. 
ADEV is the most common measure of time-domain sta- 
bility, but other related types of measures such as the 
modified Allan deviation (MDEV), Hadamard deviation 
etc are also very useful for analyzing and identifying noise 
(e.g. the MDEV in particular can distinguish between 
white and flicker phase noise). Using the Allan deviation 
it is relatively easy to identify noise processes that hap- 
pen over a long timescale. Just as in the case of the PSD 
we categorize the noise depending on the slope of the 
curved when plotted in a log-log plot. Table U shows a 
summary. Note that we have included frequency drift in 
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FIG. 5. (colour online) Allan deviation plot for the dielectric 
resonator, with(diamonds) and without (points) temperature 
stabilisation. 



table |U although stricly speaking not a noise, drift due 
to e.g. temperature fluctuations are often seen in real 
measurements. 

In order to demonstrate how ADEV can be used to 
analyze the intrinsic noise of resonators we first show the 
experimental data from the DRO in figure [5J White fre- 
quency noise (slope -0.5) dominates as expected for short 
times until linear drift (slope +1) takes over and the devi- 
ation increases rapidly. This large drift in the resonance 
frequency (about 10-20 Hz/second)is due to a changing 
temperature and agrees with the specified temperature 
stability of a few ppm/K. By employing some rudimen- 
tary temperature stabilization we can shift the onset of 
the drift further to the right on the plot. According to 
this plot the best one can hope for is a short term stability 
(and read-out accuracy 5f) of about 0.2 Hz; which agrees 
with the white noise level measured using an FFT ana- 
lyzer as well as the direct read-out using a counter. This 
corresponds to a line splitting factor Av/Sf of 5 • 10 6 . 
Figure [6] shows the ADEV -calculated using the same 
dataset as was used for the PSD- for the lumped element 
resonator measured at 30 mK. It is quite clear that the 
behaviour is very different from that of the dielectric res- 
onator. Flicker frequency noise (slope 0) dominates at 
both long and short timescales, but the dominant fea- 
ture is the sudden rise at timescales between 1 ms to 



TABLE I. Noise processes and the resulting slope in the root 
PSD and Allen deviation(ADEV) plots 



Noise type root PSD slope ADEV slope 

White phase 1 -1 

Flicker Phase 0.5 -1 

White frequency -0.5 

Flicker frequency -0.5 

Random walk -1 0.5 

Linear drift - 1 
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FIG. 6. Allan deviation plot for a 7.5 GHz lumped element 
resonator measured at 30mK. The power in the resonator is 
approximately -100 dBm. 

Is. We attribute this to random walk noise (slope 0.5). 
It is interesting to compare this with the PSD shown in 
fig. |U One obvious difference is that we can now dis- 
tiguish the types of noise much more easily. The random 
walk noise -which is the most prominent feature in fig. [5} 
is quite difficult to distinguish from other types of noise 
in the PSD and the ADEV also makes it easier to see over 
which timescale a particular noise process dominates. 

Note that the line splitting factor at long times for 
the LE resonator is about 1500; much worse than for the 
dieletric resonator. This is simply a consequence of the 
fact that the lumped element resonator is -despite having 
a much higher Q- much worse at stabilizing the Pound 
loop. This in turn can be attributed to the fact that the 
LE resonator is subject to noise from its environment: the 
flicker frequency noise is most likely due to the presence 
of two-level fluctuators^; the random walk noise might 
be due to random changes in the dielectric environment 
(due to hopping) or the movement of trapped vortices in 
the superconducting film. 

V. CONCLUSION 

By bringing well-known measurement techniques and 
analysis methods from time and frequency metrology to 
bear on the problem of the nature of noise in supercon- 
ducting resonators one can gain a lot of information not 
readily obtainable using more widely known techniques. 
By implementing a Pound loop and analyzing the data 
using the Allan deviation we are able to identify noise 
processes that are hard to distinguish by conventional 
means. Moreover, the fact that the loop directly mea- 
sures the resonance frequency makes the analysis more 
straigtforward. Finally, we would also like to mention 
that the ability to "track" the value of v$ can also lead 
to significant improvement in measurement throughput. 
A good example is measurements of v$ vs. temperature 
which can e.g. be used to determine the intrinsic loss tan- 



gent due to two-level fluctuators in the sample. This can 
be a time-consuming process if one needs to "step" the 
temperature and wait for it to stabilize before acquiring 
the scattering parameters (S21 and/or Sll) using a VNA. 
However, by using a Pound- loop one is only limited by 
the thermal time-constant between the sample and the 
temperature sensor, allowing for the temperature to be 
ramped relatively quickly. 
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